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ABSTRACT 

We report on the first results from a redshift survey of a flux-limited sample of 
X-ray clusters selected serendipitously from the ROSAT PSPC data archive. We 
spectroscopically confirm 15 clusters in the range 0.3 < z < 0.7, to a flux limit of 
~ 3.9 X 10“^^ ergs“^ cm“^, over a survey area of 17.2 deg^. The surface density of 
clusters in our survey is 2 . 0 lg '3 deg in good agreement with the number density 
of cluster candidates detected using algorithms designed to search for very extended 
sources. The number of clusters detected between 0.3 < z < 0.7 is consistent with 
a prediction based on a simple extrapolation of the local X-ray cluster luminosity 
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function, which indicates that over this redshift range no significant evolution in the 
cluster population has taken place. These results are in conflict with recent claims that 
the number density of X-ray clusters found in deep ROSAT PSPC pointings evolves 
rapidly beyond ^ = 0.3. 


Subject headings: galaxies: clusters: general — galaxies: evolution — X-rays: galaxies 
— X-rays: general 


Introduction 


Galaxy clusters are the largest gravitationally bound structures in the universe and can be 
observed out to high redshift, providing a rich source of information on the history of structure 
formation. The form of cluster evolution with look-back time is a very sensitive probe of 
hierarchical structure formation models. Assuming a CDM-like power spectrum and a self-similar 


model for structure evolution, Kaiser (1986) predicted that the comoving number density of 
clusters should increase with redshift, i.e. positive evolution. More sophisticated analytic and 
N-body models, which include gravitational instability theory and hydrodynamics, predict a 
decrease in the number density of clusters with redshift (negative evolution) for only the very 
brightest clusters (X-ray luminosities, 5 x 10^^ergs” 

luminosities (e.g. Kaiser 1991j; Cen &: Ostriker 1994). 


^44 c 1 ^ little evolution at lower 


Until recently the only available distant cluster catalogues were optically selected, showing 
little evidence for evolution in the cluster space density out to z ~ 0.5 (e.g. punn, Hoessel &: Oke| 


1986; Postman et al. 1996). However, one of the most intriguing results to emerge in recent years 


is the lack of high redshift X-ray luminous clusters in the Einstein Extended Medium Sensitivity 
Survey (EMSS, pioia et al. 1990| ; Henry et al. 1992) and other X-ray surveys (pdge et al. 1990|) . 
These studies concluded that there has been significant negative evolution over short look-back 
times in the bright end of the cluster X-ray luminosity function (XLF). The evolution found 
by Edge et al. (1990) is likely to be the result of unfortunate sampling in the redshift range 
0.1 < z < 0.15 (Ebeling et al. 1995), while the evidence for negative evolution reported by Henry 


et al. (1992) is critically examined by Nichol et al. (1997) using ROSAT pointed observations of 
EMSS clusters. 

We are now able to test the negative evolution scenario using the publicly available data 
archive of the ROSAT Position Sensitive Proportional Counter (PSPC). The greater spatial 
resolution and sensitivity, combined with the lower background, of the ROSAT PSPC compared 
to the EINSTEIN IPC means that cluster samples can be compiled to a flux limit an order 
of magnitude fainter than that of the EMSS. We have carried out a survey for serendipitously 
detected X-ray clusters from the PSPC data archive in order to study the cluster population at 
high redshift (Burke et al. 199"^). This project is referred to as the Serendipitous High-redshift 
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Archival ROSAT Cluster survey (SHARC), although our selection technique is completely different 
to that described in Nichol et al. (1997)| who use a source detection algorithm based on a wavelet 
transform to re-examine the EMSS cluster sample. There are several groups currently using the 


PSPC database to compile distant cluster surveys (RIXOS, Castander et al. 1995; RDCS, Rosati 


|et al. 19^ ; WARPS, pcharf et al. 199^ ). The purpose of this letter is to report on the number 
of spectroscopically confirmed high redshift clusters in our survey and to compare our results 
with those of other groups. In particular we make a direct comparison with RIXOS, who have 
claimed that the dearth of high redshift clusters in their sample strengthens the claims of negative 
evolution (Castander et al. 1995). In this letter we use Hq = 50kms“^ Mpc“^ and qo = 0.5. 


2. The Cluster Sample 


The ROSAT PSPC provides wide field imaging capabilities with moderate spatial and spectral 
resolution ( Triimper 19831) . On axis the PSPC has a PSF with a FWHM of 25" at 1 keV and 
an energy resolution AE/E ~ 0.4 over the range 0.1 — 2.4 keV. Although the instrument has a 
1° radius field of view, the off axis PSF degradation and the vignetting caused by the window 
support structure generally limit analysis to within a radius of 18k The X-ray emission from rich 
clusters of galaxies is spatially extended, with a surface brightness profile well fit by a King profile; 


S,,(r)(x(l + (r/r,)2)V2-3/l^ 


( 1 ) 


with a core radius (rc) of 200 - 500 kpc and /3 ~ 2/3 ( pones fc: Forman 1984 ). Assuming no 
evolution in rc and /3, clusters can be resolved out to z ~ 1.0 by the PSPC. Ideally we would like 
to optically identify all the X-ray sources in our survey, regardless of extent. However, AGN are 

(peorgantopoulos 


the primary X-ray source population at flux limits of 10 — 10 erg s ^ cm ^ 


et al. lO"^ ), outnumbering clusters by 10 to 1. Since X-ray emission from AGN is unresolved by 
the PSPC, we can significantly reduce the contamination these sources introduce into our survey 
by restricting our attention to extended sources only. Clusters have typical ICM temperatures of 2 
- 10 keV ( David et al. 199^ ), so we can reduce the contamination from soft sources and Calactic 
emission by restricting our analysis to the 0.5 — 2.0 keV pass band. 

Our survey consists of 66 deep ROSAT pointings which satisfy the following criteria: exposure 
times greater than 10 ks, |6| > 20° and 6 < 20°. The Galactic neutral hydrogen column density 
for these fields varies between 2 and 7 x 10^*^ cm“^ which, together with the average exposure 
time of 18 ks, results in an average flux limit ~ 3.9 x 10“^^ergs“^ cm“^. Due to the degradation 
of both the PSF and the detector effective area with off-axis position, we limit our survey to the 
central 18' of each field, which results in a total survey area of ~ 17.2 deg^, after masking out 
the original targets of the observations. Briefly, the analysis begins with screening out data taken 
during periods of high aspect error or particle background. We compute a global estimate of the 
background and use a sliding box technique, with a detect cell size of 30", coupled with the Cash 
statistic ( pash 1979 ) to search for sources. These sources are tested for extent by comparing the 
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photon distribution to the PSF, modelling both the positional and spectral dependence of the 
PSF. Those sources which show significant (> 3cr) angular extent are flagged as cluster candidates. 
The degree of completeness in our survey depends heavily on the extent of the cluster X-ray 
emission and the sensitivity of the extended source algorithm (e.g. see Scharf et al. 199^ . Our 
detection method is designed to pick-up moderately extended sources of size < 1' which is a good 
match to the expected angular size of clusters in the range 0.3 < z < 1. At this stage we do 
not know to what extent our completeness is compromised by ignoring intrinsically unresolved 
sources. However, both the surveys of Georgantopoulos et al. (1996) and RIXOS, which have 
optically identified all point sources in deep ROSAT PSPC fields, detect very few clusters with 
point-like X-ray emission (Georgantopoulos et al. 1996; Gastander 1996a, private communication). 
In addition, we are in the process of carrying out simulations to estimate the completeness and the 
sensitivity of our selection method to sources of different fluxes with sizes ranging from point-like 
to very extended (> 1'). 


We have carried out an optical identification programme of the extended sources, obtaining 
both imaging and multi-object spectroscopy using the European Southern Observatory 3.6m 
and the Anglo-Australian 3.9m telescopes. Typically we have redshifts for 5 — 10 galaxies per 
cluster. Our final sample consists of clusters which satisfy the joint criteria of both extended X-ray 
emission and multiple redshifts, which ensures minimal contamination from AGN. Full details of 
the survey, including the cluster identification procedure, redshift reduction and an estimate of 
the survey completeness from simulations will be given in a forthcoming paper (Burke et ah, in 
preparation). As a result of the optical observations we have identified 35 clusters ranging from 
nearby poor clusters to high redshift rich clusters. Figure || shows our four most distant clusters 
with X-ray contours overlaid on CCD R band images. 


3. Comparison with other ROSAT Surveys 


We have spectroscopically conhrmed a total of 35 clusters over an area of 17.2 deg^, giving a 
cluster surface density of 2.0j'jQ3 deg(la Poisson errors taken from Gehrels 1986| ). Two surveys 
designed to detect very extended sources in ROSAT PSPC pointings have recently published 
surface densities of extended sources found in the 0.5 — 2.0 keV band: RDCS use a wavelet 
transform algorithm and report a surface density of ~ 2.5 deg“^ at a flux limit of 3 x 10“^^ 
ergg-icm“^; WARPS use a Voronoi tesselation and percolation algorithm and detect 2.1 - 3.1 
deg“^ above an intrinsic flux limit of 4.5 x 10“^^ ergs“^cm“^. Although neither sample is 
fully identified ( Rosati et al. lO"^ report spectroscopic confirmation for a subset of their sample 
whereas Scharf et al. 1996 leave discussing their optical follow up to a future paper), their results 
are similar to our surface density of spectroscopically confirmed clusters. 


In Table we summarise our survey statistics along with those of RIXOS and show the 
corresponding redshift distributions in Figure Despite the similarity of the surveys in terms of 
limiting X-ray flux and total area surveyed, it is immediately apparent that we detect many more 
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high redshift clusters than RIXOS; beyond z = 0.3 we find 15 clusters compared to their 4. Using 
Poisson statistics, the probability that this difference could arise by chance is 2.1 x 10“^ and is 
therefore ruled out at a significance of 3.5 a. The lack of high redshift clusters in the RIXOS 
sample could be due to their detection method, as suggested by pcharf et al. 1996 , or due to 
optical misclassification. 


4. Evidence For Evolution 


We examine the redshift histogram in Figure ^ for evidence of evolution by predicting the 
expected number of clusters in the interval 0.2 < z < 0.7 using the low redshift cluster XLF of 


De Grand! (1996) . For clusters at lower redshifts, the results become unreliable due to the small 
volume sampled and the uncertainties in the XLF at low luminosities. A similar calculation was 


carried out by Castander et al. (1995) for the RIXOS survey using the cluster XLF of Edge et al. 


(1990). The cluster XLF of pe Grand! (1996)| is well matched to our dataset since it is a purely 
X-ray selected cluster sample, identified from the ROSAT All Sky Survey. The differential XLF, 
n(L), is modelled using a Schechter function. 


n(L 44 ) = Aexp (—L 44 /L*) L 44 ", 


( 2 ) 


where L 44 is the X-ray luminosity in units of 10 


i44ergs a = 1.32]']| 


-1-0.21 

0.23 


, T* = 2.631°!^ lO^^ergs 


.,-1 


and A = 4.29 ](q3® X 10 ^Mpc ^(10^'^ergs ^)" 


Using this model we can estimate the number of high redshift clusters we expect to detect 
by integrating equation ^ above the minimum detectable luminosity for each redshift interval. 
The exposure time, background count rate and a model of the off axis behaviour of the PSF 


( Hasinger et al. 1994 ) are used to predict the minimum detectable count rate at each redshift, 
after correcting for the flux falling outside the detection cell using equation ||, with /3 = 2/3 and 
Tc = 250 kpc. We use a 6 keV thermal bremsstrahlung model to convert count rates to fluxes, 
which are then corrected for Galactic absorption and k-corrected to the 0.5 — 2.0 keV band. Both 
our survey and the EMSS have a flux limit which varies from field to field. While this complicates 
the analysis, the VjVmax method used by Henry et al. (1992) fully corrects for this effect and our 
procedure follows the same technique. 

The results summarised in Table ^ demonstrate that in the redshift intervals 0.2 — 0.3 and 
0.3 — 0.7, the number of X-ray clusters predicted and observed differ by only 1.3 cr and 0.9 cr 
respectively. Despite the possibility of residual incompleteness our data provides no evidence for 


the negative evolution reported by [Castander et al. (1995) 


To test whether uncertainties in the local XLE can bias this analysis, we have repeated the 
calculation using the cluster XLF of Ebeling et al. (1996)| , which is not strictly X-ray selected, but 
based partly on the optical catalogues of Abell and Zwicky. The XLF is modelled as before, with 
a = 1.78T0.09, L* = (8.92T1.66) lO^'^ergs"! and A = (5.72T0.95) x lO-"^ Mpc-^ (10“^^ ergs”!) "“i. 
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over the energy range 0.1 — 2.4 keV. Converting to the 0.5 — 2.0 keV band, we predict 6.9l{ g 
clusters between 0.2 < z < 0.3 and 23.1^48 clusters between 0.3 < z < 0.7. These results are 
virtually identical to those obtained above, indicating that our predictions are insensitive to 
potential biases in the local XLF. 

One difference between our survey and RIXOS is that, in common with EMSS, our flux 
limit varies with the exposure time and Galactic neutral hydrogen column density of each 
observation, with an rms variation of ~ 20%, whereas RIXOS impose a constant flux limit of 
3 X 10“^^ ergs“^ cm“^ across all fields. Therefore one possibility is that the ‘missing’ clusters lie 
in our deepest fields, with fluxes below the RIXOS flux limit. However, if we replace the individual 
flux limits of each field with the average value quoted in Table [T|, the predicted number of clusters 
remains essentially unchanged, lO.T^gg compared with I9.4^|'g. Hence the deeper fields in our 
survey do not unduly bias our results. 


5. Conclusions 


We have carried out a redshift survey of extended sources selected from the ROSAT PSPC 
pointed data archive over an area of 17.2 deg^ and to a depth of 3.9 x 10“^^ ergs“^ cm“^ in 
the energy range 0.5 — 2.0 keV. Spectroscopic follow-up observations have produced a sample of 
distant X-ray selected galaxy clusters. The surface density of our spectroscopically confirmed 
cluster sample is 2 . 01 q ;3 deg ™ good agreement with the number density of cluster candidates 
detected using algorithms designed to search for very extended sources (RDCS and WARPS). 

In addition, the redshift histogram of our clusters is consistent with a simple extrapolation of 
the local luminosity function of X-ray clusters from De Grand! (1996) . This indicates that no 
substantial negative evolution has taken place in the cluster population over the redshift range 
0.3 < z < 0.7. 


After this manuscript was accepted we spectroscopically confirmed a further cluster 
at z = 0.672 from service observations carried out on the 4.2m WHT. This completes the 
identification of extended sources and brings the number of clusters in our survey which are 
beyond z = 0.3 to 16, which marginally increases the discrepancy between our sample and that of 
RIXOS. 


This research has made use of data obtained from the Leicester Database and Archive Service 
at the Department of Physics and Astronomy, Leicester University, UK. We thank Sabrina De 
Grand! for providing us with the errors on the ESO KP cluster XLE. Mel Ulmer and Brad Holden 
are acknowledged for their support of the SHARG survey. GAG and DJB acknowledge PPARG 
for an Advanced Eellowship and Postgraduate Studentship respectively. AKR acknowledges 
support from NASA ADP grant NAG5-2432. We thank the staff at both the European Southern 
Observatory and the Anglo-Australian Observatory for use of these facilities. 
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Fig. 1.— X-ray contours, for the 0.5 — 2.0 keV band overlaid on CCD R band exposures of our 
4 most distant clusters. The images are 10 minute exposures taken at the ESO 3.6m. Axis units 
are in arcminutes and the X-ray images have been lightly smoothed and begin at 3a above the 
background with a step size of 3a. A is at z = 0.576, B at z = 0.550, C at z = 0.505 and D at 
z = 0.484. At these redshifts 1 arcminute subtends ~ 440 kpc. 


Fig. 2.— Comparison of our complete cluster redshift distribution (solid line) with that of the 
RIXOS survey (dotted line). The RIXOS data are taken from Castander et al. (1995) for z > 0.2 
and Castander (1996b) for z < 0.2. 
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Survey O Flux limit N^. Nc 

(deg^) (10“^'^ ergs“^ cm“^) (0.2 < 2 ; < 0.3) (0.3 < 2 ; < 0.7) 

This Work IF2 T9 15^^ (lO.dt^;^) 

RIXOS 14.9 3.0 Otlg (8) (25) 


Table 1: Comparison between this work and RIXOS for all 2 ; > 0.2 clusters. The errors on the 
number of clusters found are Poisson. The numbers in brackets for our survey are the predictions 
for the number of clusters observed based on a model of our selection function and the XLF of De 
Grandi (1996), along with the Icr errors from the simulations. The values for RIXOS have been 
taken from Castander et al. (1995), which uses the XLF of Edge et al. (1990). 
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